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ABSTRACT 


A model  was  developed  to  combine  radiosonde  observations  with  ultra* 
sensitive  well-calibrated  radar  measurements  of  clear  air  turbulence  to 
provide  estimates  of  two  eddy  parameters,  the  eddy  dissipation  rate  and 
the  eddy  heat  diffusion  coefficient.  The  eddy  parameters  were  obtained 
as  volume  averages  over  1 km  thick  layers  and  time  averages  over  several 
liours.  The  model  was  tested  using  simultaneously  obtained  radar  measure- 
ments and  pilot  reports  of  turbulence  intensity  from  a U-2  aircraft.  The 
model  was  also  shown  to  be  consistent  with  other  reported  radar  and  air- 
craft measurements. 

Simultaneous  radar  and  radiosonde  measurements  were  obtained  on  24 

different  days  spanning  all  seasons.  The  median  eddy  heat  diffusion 

coefficient  for  all  the  observations  and  the  13  to  16  km  height  region 
2 2 

was  0.02  m /sec  (200  cm  /s) . This  value  did  not  change  from  season  to 
season. 
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1.1  Objectives 

The  objective  of  this  contract  is  the  estination  of  turbulent  eddy 
diffusion  coefficients  for  the  lower  stratosphere  using  a collection  of 
simultaneous  radar  and  radiosonde  observations.  An  important  aspect  of 
the  study  reported  herein  is  the  use  of  a unique  set  of  radar  data 
obtained  by  MIT  Lincoln  Laboratory  at  the  Millstone  Hill  Radar  Facility 
during  the  1968-1970  time  period. 

1.2  Summary  of  Results 

The  simultaneously  acquired  radar  and  radiosonde  data  were  processed 
to  provide  estimates  of  the  turbulent  eddy  dissipation  rate  and  eddy 
heat  diffusion  coefficient.  Median  profiles  of  these  eddy  parameters 
are  given  in  Table  I.  The  profiles  represent  the  median  values  of  the 
medians  for  each  observation  set.  Thirty-two  sets  of  observations  were 
obtained  on  24  different  days:  3 days  between  11  January  and  10  February 
1968;  4 between  28  May  and  27  June  1968;  4 between  29  July  and  2 August 
1968;  7 between  20  July  and  31  July  1970;  and  6 between  26  October  and 
12  November  1970. 

The  data  processing  started  with  profiles  of  the  refractive  index 
2 

structure  parameter  , produced  by  the  well  calibrated  Millstone  Hill 

L-band  radar  (Crane,  1970;  Crane,  1974).  The  C ^ profiles  were  trans- 

n 2 

formed  into  profiles  of  the  temperature  structure  parameter,  C_  , using 

2 * 

the  radiosonde  data.  A model  was  developed  to  use  the  C.^.  estimates  and 
radiosonde  data  to  provide  the  required  eddy  parameter  estimates  plus 
estimates  of  the  volume  filled  by  turbulence.  The  model  was  verified 
using  pilot  reports  of  turbulence  obtained  during  U-2  aircraft  flights 
made  in  conjunction  with  5 of  the  observation  sets. 

The  radar  deprived  eddy  parameters  represent  one  cf  the  most  exten- 
sive data  sets  available  to  date  on  the  intensity  of  small  scale  lower 
stratospheric  turbulence  at  a single  geographical  location.  The  eddy 
heat  diffusion  coefficient  estimates  presented  in  Table  I for  the  13-16  km 
height  interval  are  identical  to  the  mean  value  computed  by  Lilly  et  al 
(1974)  based  upon  Project  HICAT  data  for  the  13.7  - 15.2  km  height  interval 

7 


c 

o 


<0 

U) 

tA 

tn 

lO 

1/) 

D. 

1 

1 

1 

1 

1 

t/) 

O 

O 

o 

o 

o 

</>  CJ 

t/t  ♦JCM 

(/) 

QC 

(0  £ 

Q OC  w 

K 

K 

X 

X 

K 

UJ 

H 

UJ 

dy 

»0 

ro 

hO 

<M 

CQ 

-t 

H 


< 

o. 

>• 

o 

o 

U4 

z 

< 

o 

UJ 


•o 

LU 


c 

o 


3 

4-> 

C 

U4 

O 

rj 

rj 

CM 

CM 

CM 

•H 

1 

1 

1 

1 

1 

o 

o <N» 

O 

c 

O 

O 

o 

4-» 

6 

<4-*  W 

(0 

U-4 

X 

X 

X 

X 

X 

OJ 

z 

S|J= 

u 

CM 

CM 

CM 

ro 

•o 

UJ 


w 

JZ 

> 

rO 

in 

o 

r^ 

00 

£ 

«-H 

• M 

0) 

1 

1 

1 

1 

1 

4J 

CM 

x> 

m 

vO 

X 

c 

and  similar  underlying  terrain.  The  12-13  km  values  in  Table  I are  also 
consistant  with  the  estimates  of  mean  eddy  parameters  reported  by  Cadet 
(1977)  although  his  data  are  for  a single  set  of  observations  over  moun- 
tainous terrain.  The  eddy  diffusion  coefficients  are  an  order  of  magni- 
tude smaller  than  the  estimates  provided  by  Rosenberg  and  Oewan  (1974). 
Their  estimation  procedure  employed  however,  a free  parameter,  the  time 
stop,  which  must  be  adjusted  to  match  the  eddy  dissipation  rate. 

The  results  presented  herein  are  consistant  with  the  available  data 
on  transport  by  small  scale  turbulence  in  the  lower  stratosphere.  They 
are  onc-to-two  orders  of  magnitude  smaller  than  the  coefficients  used  to 
adjust  the  output  of  one-dimensional  modeling  analyses  to  match  observa- 
tions of  trace  constituent  distributions.  This  large  discrepancy  arises 
because  the  onc-dimensional  models  are  generally  formulated  to  use  the 
turbulent  diffusion  model  to  account  for  larger  scale  non-diffusive  trans- 
port (sec  for  example,  Reiter,  1975). 
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2.1  The  Roles  of  Stratospheric  Turbulence 

Initially,  interest  in  clear  air  turbulence  (CAT)  in  the  lower 
stratosphere  and  upper  troposphere  arose  because  of  its  potential  hazard 
to  aircraft.  Project  HICAT  (Crooks  et  al,  1967)  was  conducted  by  the 
U.S.  Air  I'orce  to  establ ish  structural  design  criteria  for  aircraft.  The 
single  U-2  aircraft  used  for  project  HICAT  was  equipped  to  obtain  true 
gust  velocity  data  on  scale  sizes  ranging  from  20  meters  (m)  to  20  kilo- 
meters (km).  The  HICAT  U-2  flew  over  6.6  x 10^  km  sampling  turbulence 
in  the  13  to  20  km  height  range.  These  data  are  the  most  extensive 
available  for  the  study  of  CAT  in  the  lower  stratosphere.  Lilly  et  al 
(1974)  analyzed  these  data  to  provide  eddy  heat  diffusion  coefficient 
est irnates  for  several  height  regions  and  different  underlying  terrain 
conditionh.  The  HICAT  instrumentation  was  designed  to  obtain  true  gust 
velocities  - velocity  variance  data  on  the  100-300  m scale  size  of  impor- 
tance to  aircraft  design.  The  instrumentation  was  not  designed  however, 
to  obtain  the  data  at  the  much  smaller  scale  sizes  within  the  inertial 
subrange  of  the  turbulence  spectrum  required  for  the  direct  estimation 
of  eddy  parameters. 

Currently,  interest  in  CAT  in  the  lower  stratosphere  is  centered 
on  the  estimation  of  vertical  eddy  diffusion  coefficients  for  use  in 
modeling  the  effects  of  stratospheric  pollutants.  The  stable  region 
between  1 < and  16  km  height  acts  as  a barrier  to  the  vertical  transport 
of  f luorocarborns  from  the  lower  atmosphere  into  the  stratosphere  where 
they  participate  in  catalytic  reactions  to  reduce  ozone  concentrations 
(Molina  and  Rowland,  1974).  The  lower  stratosphere  also  acts  as  a barrier 
to  the  vertical  transport  of  long  lived  foreign  gasses.  This  latter 
problem  has  been  actively  studied  under  the  sponsorship  of  the  U.S.  Dept, 
of  Transport  to  assess  the  effects  of  supersonic  transport  (SST)  air- 
craft (see  the  CIAP  monographs,  1975). 

Models  generated  for  the  analysis  of  the  effects  of  possible  stra- 
tospheric pollutants  depend  critically  on  the  values  assumed  for  the 
eddy  diffusion  coefficient  . These  models  must  parameterize  the  role  of 
subgrid  scale  transport  processes.  One-dimensional  models  must  combine 
the  effects  of  small  scale  turbulence,  larger  scale  transport  by  convective 
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elements,  transport  by  larger  scale  folds  or  breaks  in  the  tropopause 
(Danielson,  1968)  and  perhaps  other  larger  scale  processes.  Of  these, 
only  transport  by  small  scale  turbulence  is  truly  diffusive  and  amenable 
to  analysis  using  the  gradient  diffusion  equation.  The  selection  of  a 
diffusion  coefficient  for  a one-dimensional  model  is  usually  an  exercise 
in  curve  fitting  using  observations  of  a tracer  (see  for  example,  Reiter 
and  Bauer,  197S) . The  same  may  also  be  true  of  three-dimensional  models 
if  the  larger  scale  processes  are  not  included  explicitly. 

A study  of  small  scale  turbulence  is  important  to  the  understanding 
of  the  possible  deficiencies  of  a model  and  to  provide  estimates  of  the 
minimum  values  to  be  expected  if  the  other  larger  scale  processes  are 
not  operative.  Unfortunately,  the  intensity  of  turbulence  is  not  a para- 
meter of  the  medium  but  a parameter  of  the  larger  scale  atmospheric 
motions.  The  eddy  diffusion  coefficients  are  not  constant  but  depend 
upon  the  larger  scale  processes.  The  studies  of  small  scale  turbulence 
therefore  must  be  morphological  - they  must  establish  both  the  magnitude 
of  the  sm.*!!!  scale  transport  and  its  variation  and  their  dependences  on 
the  large  scale  flow. 
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2.2  Radar  as  a Tool  for  the  Measurement  of  Stratospheric  Turbulence 

Studies  of  small  scale  turbulence  in  the  lower  stratosphere  can  be 
performed  using  insitu  observations  from  aircraft  or  balloon  borne  plat- 
forms or  remotely  using  monostatic  or  bistatic  radar  systems.  The  air- 
craft measurements  (Crooks  et  al,  1967)  or  balloon  measurements  (Cadet, 
1977)  currently  available,  however,  do  not  sense  fluctuations  at  small 
enough  scales  to  provide  direct  observations  of  the  eddy  parameters. 
Monostatic  radar  systems  can  directly  observe  the  intensity  of  the  tur- 
bulent refractive  index  fluctuations  at  scale  sizes  equal  to  one-half  the 
radar  wavelength.  These  observations  must  also  be  interpreted,  however, 
to  provide  eddy  parameter  estimates. 

Radar  observations  of  CAT  have  been  available  for  over  a decade. 
Hardy  and  Katz  (1969)  reviewed  the  observations  of  CAT  and  showed  that 
high  powered  radar  systems  were  useful  for  morphological  studies  of  the 
occurrence  of  thin  turbulent  layers  in  the  troposphere.  Occasionally, 
radar  detections  of  the  tropopause  were  possible  (Atlas  et  al , 1966). 
Using  the  high  powered  Millstone  Hill  L-band  radar.  Crane  (1970)  reported 
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the  regular  detection  of  thin  turbulent  layers  at  the  tropopause  and  at 
heights  up  through  20  km.  On  one  occasion  a layer  was  detected  at  a 
height  of  25  km. 

In  the  early  radar  observations,  the  main  interest  was  to  demon- 
strate that  the  echo  structure  observed  by  the  radar  was  indeed  CAT  and 
was  correlated  with  CAT  as  observed  by  aircraft.  As  the  association 
between  CAT  and  the  observed  radar  echoes  became  better  established,  the 
interest  changed  to  obtaining  a better  picture  of  the  production  mecha- 
nism. The  observations  of  the  marine  inversion  layer  reported  by  Atlas 
et  al  (19"0)  and  by  Metcalf  and  Atlas  (1973)  using  the  extremely  high 
resolution  FM-CW  radar  and  observations  of  large-amplitude  Kelvin-Helm- 
holtz  billows  reported  by  Browning  (1971)  showed  that  the  thin  turbulent 
layers  occur  in  stably  stratified  shear  flow  when  the  gradient  Richardson 
number  is  less  than  0.25.  The  structure  of  the  resultant  turbulence 
traced  the  classical  pattern  for  the  initial  stages  of  roll-up  by  Kelvin- 
llclmholtz  instability.  The  association  between  the  thin  layer  radar 
echoes  and  turbulence  was  complete;  the  mechanism  for  CAT  in  stably  stra- 
tified shear  flow  revealed. 

Currently,  radar  observations  of  lower  stratospheric  CAT  are  used 
to  obtain  wind  profiles  (Green  et  al,  1975;  Balsley  and  Farley,  1976). 

The  turbulent  fluctuations  of  refractive  index  are  used  as  tracers  for 
the  measurement  of  the  mean  wind  profile.  Radar  data  have  not  been  used 
to  provide  eddy  parameter  estimates  because  most  of  the  radars  are  not 
adequately  calibrated  and  an  adequate  analysis  model  was  not  available. 

Crane  used  the  high  powered  Millstone  Hill  L-band  radar  to  obtain  well 
calibrated  radar  measurements  of  CAT  in  the  upper  troposplere  and  lower 
stratosphere  as  an  element  of  a larger  propagation  measurement  program 
conducted  for  NASA  (Crane,  1973)  and  for  the  U.S.  Air  Force  (Crane,  1970). 

The  object  of  this  contract  for  the  Air  Force  Geophysics  Laboratory  is 
to  develop  a model  for  the  interpretation  of  the  radar  data  and  to  repro- 
cess the  data  acquired  between  1968  and  1970  to  obtain  eddy  diffusion 
coefficient  estimates. 
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3.  RADAR  MEASUREMENTS  OF  STRATOSPHERIC  TURBULENCE 


3.1  Radar  Estimation  of  the  Refractive  Index  Structure  Function 

Fluctuations  of  the  index  of  refraction  within  a turbulent  region 
can  produce  scattering  of  electromagnetic  waves  propagating  through  the 
region.  Sensitive  radars  can  detect  the  scattering  to  remotely  sense 
the  intensity  of  the  refractive  index  fluctuations.  Tlie  theory  of  scat- 
tering of  electromagnetic  waves  by  a turbulent  medium  is  well  established 
(Tatarskii,  1967).  The  bistatic  radar  cross  section  per  unit  volume  of  a 
turbulent  medium  is  related  to  the  three-dimensional  (3-d)  power  spectral 
density  of  refractive  index  fluctuations,  by  (Crane,  1968;  Crane,  1973): 


B = Sit^k  m<t'^(J<-kp) 


(1) 


where  P 
k 
m 


n 


= bistatic  scattering  cross  section  per  unit  volume 
= 2it/X  = wavenumber;  X = wavelength 
= polarization  mismatch  factor 

= 3-d  power  spectral  density  of  refractive  index  fluctuations 

= kx  where  x is  a unit  vector  from  the  transmitter  to  the 
scatterer 

= unit  vector  from  the  scatterer  to  the  receiver. 


Bistatic  refers  to  a radar  system  with  the  transmitting  antenna  physi- 
cally separated  from  the  receiving  antenna.  When  the  same  antenna  is 
used  for  both  transmission  and  reception,  the  system  is  called  a mono- 
static radar  or,  more  commonly,  a radar. 

For  clear  air  turbulence  sampled  at  a wavenumber  Ijk-kp]  in  the 
inertial  subrange  (Tatarskii,  1967): 


4>  (s)  = 4>  (s)  = 0.033  C ^ 
n 'C  n ' n 


(2) 


2 » « » 
where  C^  is  the  refractive  index  structure  parameter  and  s = = k(x-p) 

The  radar  observable,  the  scattering  cross  section  per  unit  volume 

therefore  is  related  to  the  refractive  index  structure  parameter  by: 


p = 8TT2k‘^4>^(2)<)  = 


(3) 
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where,  for  the  inonostatic  radar  case,  p = -x,  and  m = 1 . 

A radar  observes  the  total  scattering  cross  section  obtained  by  inte- 
grating the  scattering  cross  section  per  unit  volume  over  the  resolution 
volume  of  the  radar.  The  radar  equation  for  a distributed  target  such  as 
turbulence  is  given  by 


I’t  (4it)^ 


/gl(;|c)g2(lc)e(;iC)H(jc) 


where  is  the  measured  received  power 

P^.  is  the  transmitted  power 

8.1,112  are  receiver  and  transmitter  line  losses 

gl,g2  are  normalized  radiation  patterns  for  the  receiver  and 

transmitter  antennas  (directivity) 

Gi,G2  are  the  receiver  and  transmitter  antenna  gains 

and  H(x)  is  the  range  weighting  function. 


Assuming  that  the  turbulence  uniformly  fills  the  resolution  volume,  equa- 
tion (4)  may  be  reduced  to 


(4ir)3x'* 


e(x)v(;() 


where  V(^)  = / gig2H  is  the  resolution  volume  for  a radar  return  at 


Using  the  equations  (3)  and  (5),  the  radar  observables  P and  X may  be 

2 r 'V, 

converted  into  a measurement  of  C : 

n 


2 2 
C = K*P  'X 
n r 


where  X is  the  range  to  target  and  K is  the  radar  calibration  constant. 

3.2  A Model  to  Estimate  the  Intensity  of  Turbulence  Using  Radar  Data 

2 

3.2.1  The  Temperature  Structure  Parameter,  C.j, 

2 

The  radar  system  provides  an  estimate  of  averaged  over  the 
radar  resolution  volume.  The  estimate  is  derived  from  an  observation 
made  at  a spatial  wavenumber  of  2k  or  a scale  size  etiual  to  one  half  the 
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radar  wavelength,  a scale  size  that  generally  lies  in  the  inertial  sub- 
range for  the  radars  used  for  turbulence  analysis.  The  parameters  of 

interest  are  the  eddy  parameters,  the  eddy  dissipation  rate  and  eddy 

2 

heat  diffusion  coefficient,  not  C 

At  microwave  frequencies,  the  index  of  refraction  is  related  to 
the  meteorological  properties  of  the  medium  by  (Crane,  1976): 


n = 1 + 


,77.6P 


3.73  X lO^e,  ,„-6 
7 ) X 10 


(7) 


where  n is  index  of  refraction 

P is  pressure  (in  millibars,  mb) 

1 is  temperature  (in  degrees  Kelvin,  °K) 
and  e is  water  vapor  pressure  (in  mb). 


i 

t 


1 


I 

L 
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The  refractive  index  fluctuations  are  caused  by  turbulent  variations  in 
temperature  (or  density)  and  water  vapor  concentration.  At  the  heights 
of  interest  in  the  lower  stratosphere  and  upper  troposphere,  the  effects 
of  water  vapor  are  negligible.  The  relevant  equation  for  refractive  index 
variation  about  the  mean  value  reduces  to 


,77.6PA0,  ,„-6 

An  = ( T ) X 10 

T 


(8) 


where  A denotes  change  and  0 is  the  potential  temperature.  The  expression 
2 

for  therefore  may  be  simply  related  to  the  temperature  structure 
parameter,  by 


C ^ = f-X X 10^)  ^ 

^T  ^77. 6P  ^ J ^n 


(9) 


hence,  C.^,  may  be  obtained  directly  from  the  radar  observations  of  P^  and 
X when  the  pressure  and  temperature  in  the  scattering  volume  are  known. 

3.2.2  Eddy  Parameters  and  K^^ 

The  radar  may  be  used  to  obtain  an  estimate  for  the  temperature 
structure  parameter  averaged  over  the  radar  resolution  volume.  The  large 
number  of  available  radar  observations  of  CAT  show  that  the  actively  tur- 
bulent regions  are  very  thin,  and  may  in  general  be  thinner  than  the 
resolution  volume  of  the  radar  (Hardy  and  Katz,  1969;  Atlas  et  al,  1970). 
Let  the  actively  turbulent  region  occupy  a thin  layer  of  height  increment 
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d within  a radar  resolution  volume,  of  height  increment  h.  Noting  that 
the  layers  are  generally  large  in  horizontal  extent,  assume  that  the 
radar  resolution  volume  is  smaller  in  the  horizontal  than  the  actively 
turbulent  region.  Assume  also  that  the  turbulence  is  of  uniform  inten- 
sity throughout  its  active  region.  Then, 


(10) 


2 2 
where  is  the  value  measured  by  the  radar  and  C_  is  the  actual 

1 c T 

value  for  the  turbulent  region. 

Similarly,  assume  that  only  turbulence  within  the  active  region 
contributes  to  the  eddy  dissipation  rate. 

■r  _ <lEmm  (h-djemn 
">  h 


where  is  the  mean  (in  time  and  space)  eddy  dissipation  rate  for  the 

layer  of  thickness  h and  within  the  horizontal  dimensions  of  the  radar 

resolution  volume,  is  the  dissipation  rate  within  the  active  region, 

e is  the  dissipation  rate  outside  the  active  region  and  e >>  e . Then, 
mo  ° mm  mo 

^m  “ (it)  ^mm . ( ^ ^ ^ 


Similarly,  the  mean  eddy  diffusion  rate  may  be  expressed  as 


(12) 


The  region  of  the  atmosphere  of  interest,  the  lower  stratosphere  and 
upper  troposphere,  are  statically  stable.  In  such  a region,  the  eddy 
diffusion  parameters  for  momentum  and  heat  may  be  defined  by  (Monin  and 
Yaglom,  1975;  Dutton,  1969) 


^m  ” ^ (1"^) 


(13) 


and 


E.n  = 


(14) 


where  c ,c 
m q 


are  the  mean  eddy  dissipation  rate  and  eddy  heat  diffusion 
rate  values  respectively, 


K are  the  eddy  diffusion  coefficients  for  momentum  and  heat 
respectively. 


AV  is  the  mean  vector  velocity  change  (difference)  across  the 
^ layer  of  thickness  h. 
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AO  is  the  mean  potential  temperature  across  the  layer 


and 


J- 

0 

o 


* 7“ 
Z 

i = l 


CW<))> 


= WVi>(& 


A0/h 

(AJ^)Vh? 


g is  the  acceleration  due  to  gravity, 

0^  is  the  mean  potential  temperature  of  the  layer. 


<w<(i>  is  the  turbulent  heat,  vertical  velocity  covariance. 


(15) 


c p <w<()>  is  the  vertical  turbulent  heat  flux  and 
P o 

<wVj>  are  the  turbulent  velocity  covariances  with  i •=  1,2  horizontal 
directions;  is  the  vertical  momentum  transport. 


The  turbulent  kinetic  energy  balance  equation  (13)  was  derived  using 

the  standard  assumptions  of  steady  state  turbulence,  horizontal  homogeneity, 

vertical  volume  boundaries  not  within  the  actively  turbulent  region. 

The  last  two  assumptions  are  required  to  eliminate  terms  including  the 
divergence  of  the  turbulent  transport  of  turbulent  kinetic  energy  and  the 
divergence  of  the  turbulent  velocity,  pressure  covariance.  It  is  noted 
that  in  most  measurements  of  turbulence,  these  assumptions  are  often 
violated  when  the  turbulence  is  not  homogeneous  but  in  radar  measurements 
of  apparently  steady  state,  homogeneous  thin  horizontal  turbulent  layers 
thinner  than  a resolution  element,  the  assumptions  are  not.  Similar  assump- 
tions are  required  for  the  simplification  of  the  turbulent  heat  diffu- 
sion rate  equation  (14).  Both  equations  express  the  steady  state  balance 
between  larger  scale  production  and  small  scale  dissipation  assumed  for 
the  existence  of  an  inertial  subrange. 

The  turbulent  momentum  and  heat  diffusion  rate  equations  (13)  and 
(14)  do  not  require  the  existence  of  eddy  diffusion  coefficients  since 
they  were  initially  derived  using  turbulent  heat  and  momentum  fluxes. 

They  describe  volume  averaged  turbulence  production  and  dissipation. 
Similarly,  the  flux  Richardson  number  is  a ratio  of  volume  averaged 
covariances  divided  by  the  mean  wind  shear  across  the  layer.  The  average 
intensity  of  the  turbulence  within  the  volume  may  be  characterized  by 
the  eddy  parameter  the  average  vertical  transport  through  the  layer  by 
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the  eddy  heat  diffusion  coefficient  Kj^.  The  flux  Richardson  number  is 
related  to  the  gradient  Richardson  number  by 


R . ^ R.  - !i 


where  R . = 
1 


-f_  A0/h  _ _g_  A0h 
®o  (AV)2/h2  ■ Qq  (Ay)2 


is  the  gradient  Richardson 


number  and  is  the  eddy  Prandtl  number. 

The  parameters  defined  by  equations  (13)  and  (14)  are  for  a layer 
of  tliickness  h.  For  the  actively  turbulent  layer  of  thickness  d,  a new 
set  of  parameters  are  required  each  of  which  apply  only  to  the  active 


region: 


Av 

e = K (^)^(l-R-) 
mm  m d ^ 


K (^)2 


where  Ay  expresses  the  mean  vector  velocity  change  across  the  thickness 
d and  A4‘  the  potential  temperature  difference  across  the  layer.  The 
quantities  K^,  and  will  in  general  be  different  from  the  h-layer 
quantities  K^,  Kj^  and  R^.  Equations  (10)  through  (18)  are  consistent 
with  current  ideas  of  intermittent  turbulence  (Monin  and  Yaglom,  1975). 
They  do  not  provide  a complete  set  of  equations  for  the  computation  of  d 
or  any  of  the  d-layer  or  h-layer  eddy  parameters.  The  h-layer  gradient 
values  may  be  obtained  from  radiosonde  observations;  the  d-layer  values 
cannot. 

2 

An  additional  equation  relating  c , e , and  C_  can  be  obtained 
^ ^ qm’  mm  T 

from  the  classical  dimensional  analysis  of  inertial  subrange  turbulence 
(Monin  and  Yaglom,  1975) 


= a-E 
0 mm 


whcj e is  an  empirical  constant  currently  estimated  at  3.3.  Further 
prof ress  in  the  reduction  of  these  equations  requires  new  assumptions. 
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The  first  new  assumption  is  that  every  volume  in  the  atmosphere  is  on 
average  (in  time)  always  turbulent.  This  assumption  is  consistant  with 
recent  observations  such  as  those  reported  by  Goodman  and  Miller  (1977) 
which  provide  evidence  for  turbulent  mass  transport  across  a stably 
stratified  layer  with  a gradient  Richardson  number  near  24.  For  turbu- 
lent transport  across  a stable  layer,  the  time  averaged  flux  Richardson 
number  for  the  layer  must  lie  in  the  range  between  0 and  1.  If  this 
layer  is  turbulent  but  with  different  magnitudes  of  turbulence  in  sub- 
layers within  the  layer,  the  time  averaged  R^  should  be  within  this  range 
for  each  of  the  sublayers.  This  model  agrees  with  the  accepted  model  for 
turbulence  production  in  a stable  region  with  shear  flow  and  is  consis- 
tant with  available  radar  observations  and  laboratory  measurements.  A 
cascade  of  Kelvin-Helmholtz  instabilities  occur  at  random  in  time  within 
each  layer.  The  stably  stratified  shear  flow  becomes  locally  unstable 
when  the  local  gradient  Richardson  number  is  less  than  1/4.  Theoretical 
estimates  reviewed  by  Howard  and  Maslowe  (1973)  shows  that  stability  (no 
turbulence)  may  be  assumed  only  if  the  local  R^  is  greater  than  1/4  for 
all  d.  Experimental  data  reported  by  Thorpe  (1973)  showed  the  local  R^ 

< 0.2  across  the  active  region  at  the  onset  of  turbulence  and  the  local 
R.  = 0.3  at  the  cessation  of  turbulence.  A model  consistant  with  these 
observations  is  that  the  flux  Richardson  number  is  everywhere  constant 
throughout  the  layer  at  a value  between  0.2  and  0.3.  The  flux  Richardson 
number  may  be  taken  to  equal  the  gradient  Richardson  number  within  the 
actively  turbulent  region;  the  h-layer  gradient  Richardson  number,  R^  is 
given  by  the  radiosonde  derived  value  for  a thickness  h.  For  an  actively 
turbulent  layer,  the  local  R^  is  taken  as  equal  to  the  critical  value, 

R^  = 1/4. 

Using  the  concept  of  a flux  Richardson  number  defined  for  a layer. 


and 
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Si  ^/'^f 
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(aV)^  (1-Rf) 
h2  Rf 


h 


(1-Rf) 


(21) 
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For  = 1/4  = constant,  {1-R^)/R^  = 


'm  ■ ^'Sl  f f ■ <“> 

o 

where  N = Brunt  Vaisala  frequency.  The  result,  equation 

(22)  is  the  relationship  Lilly  et  al  (1974)  used  to  reduce  the  Project 
IIICAT  data  and  Cadet  (1977)  used  to  reduce  his  balloon  data.  Using 
equations  (21)  and  (11)  and  the  assumed  constancy  of  R^  over  layers  of 
thicknesses  d and  h,  R^  = R^  and 


7 -7  1 If  M g 1-R 

m"^h0  R-  ~ h ^ d e 

of  of 


or  A0  = Kj^  A41 


(23) 


and,  since  = K„  (f,=  . ^ (M,2 

A0  _ Aiji 

h ” d 

This  result  shows  that  the  potential  temperature  gradient  is  maintained 
across  the  turbulent  layer,  a result  experimentally  obtained  in  the 
laboratory  by  Long  (1973)  and  evident  in  the  aircraft  data  presented  by 
Metcalf  and  Atlas  (1973). 

By  the  definition  of  R^  and  the  assumed  equality  of  R^  and  R^  with 
the  actively  turbulent  layer  (P^  = 1 in  the  actively  turbulent  layer). 


r^i  ^ - _S_  ^ . _L 

'■  d ^ " 0 d R - 

o f 


(24) 


where  R^  is  defined  over  the  layer  of  thickness  h.  Equation  (24)  ex- 
presses the  relationship  between  the  velocity  gradients  across  the 
actively  turbulent  layer  and  the  layer  of  thickness  h.  To  close  the  full 
set  of  equations,  an  additional  telaltionship  between  A)^  and  A)^  is  required. 
Noting  the  tendency  for  stably  stratified  fluids  to  form  horizontal 


2 2 

internal  boundary  layers,  the  second  new  assumption  that  (Av)  * (A)^)  is 
made.  Using  this  assumption  and  equation  (24), 


d = 


(25) 


Equation  (25)  matches  the  equation  Long  (1973)  obtained  to  fit  his 
experimental  data  for  the  thickness  of  the  actively-turbulent  steady 
state  shear  layer. 

The  two  new  assumptions  required  to  close  the  equations  were 
argued  on  physical  grounds  and  shown  to  match  the  experimental  results 
obtained  by  Long.  Using  these  assumptions  and  equations  (10)  through 
(25),  the  final  results  for  the  analysis  of  the  radar  data  are  obtained: 


•Si  - (Se^ 


2,3/2 


1/2  ,A0r5/2  -3/2 
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(26) 
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23/2  £.5/2  ,l-Rf,3/2,h,l/2,A0,-3/2.  -3/2 

^^Te^  ^d^  “0 


(27) 


These  results  are  similar  to  those  recently  presented  by  Van  Zandt  et  al 
(1977)  and  have  definite  advantages  over  previous  analyses  (see  Hardy 
and  Katz,  1969;  Metcalf  and  Atlas,  1973)  in  that  an  outer  scale  for  tur- 
bulence does  not  have  to  be  assumed.  The  results  still  are  an  approxi- 
mation in  that  the  h-layer  contains  many  active  sublayers,  not  one  - 
each  turbulent  region  is  intermittent.  If  the  eddy  dissipation  rate 
values  for  each  sublayer  of  intermittent  turbulence  were  identical,  the 
empirical  coefficient  would  account  for  the  non-linear  averaging  of 
the  randomly  intermittent  patches  within  the  layer  of  thickness  d (or  h) 
(sec  Monin  and  Yaglom,  1975,  section  25).  Since  e can  vary  from  sub- 
layer  to  sublayer,  may  differ  from  the  assumed  value.  Noting  that 
the  larger  scale  motions  responsible  for  generating  the  turbulence  will 
not  selectively  pick  the  same  sublayers  when  h is  the  order  of  1 km,  it 
is  expected  that  the  differences  are  small  when  long  time  averages  are 
used. 
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The  final  equations  (26)  and  (27)  follow  directly  from  the  defini- 

2 

tions  of  and  and  the  inertial  subrange  expression  for  - equation 
(19)  when  h = d and  the  entire  layer  is  turbulent.  In  this  case,  the 
new  assumptions  used  to  find  d are  not  required.  The  impact  of  the  new 
assumptions  is  to  provide  a derivation  of  the  correction  factor  for 
layer  thickness,  the  factor  including  h/d  which  adjusts  the  radar  data 
because  the  thin  turbulent  layers  do  not  normally  fill  the  radar  resolu- 
tion volume.  If  a radar  with  a sufficiently  small  resolution  volume  is 
used,  such  as  the  FM-CW  radar  of  Richter  (1969),  with  a resolution  volume 
entirely  filled  by  the  active  region,  then  the  new  assumptions  are  not 
required.  In  this  case,  the  application  of  equations  (26)  and  (27), 
however,  may  be  more  difficult  to  justify  because  a small  active  region 
that  completely  fills  a small  resolution  volume  may  be  only  one  of  the 
many  intermittent  patches  within  the  larger  turbulent  region  and  the 
steady  state  and  spatially  homogeneous  assumptions  used  to  derive  the 
turbulent  energy  balance  equation  (13)  may  be  violated. 


4.  THE  RADAR  OBSERVATIONS 


4.1  The  Millstone  Hill  1,-band  Radar  System 

The  radar  system  and  its  calibration  for  use  in  observing  meteorolo- 
gical targets  has  been  described  by  Crane  (1971;  1973).  The  character- 
istics of  the  radar  as  configured  in  1968  and  1970  are  given  in  Table  II. 
Tbe  radar  used  an  84-foot  antenna  that  was  computer  controlled.  The 
antenna  gain  was  measured  using  radio  star  flux  measurements.  The  antenna 
pattern  information  required  to  compute  the  radar  resolution  volume  was 
obtained  both  for  the  transmission  and  reception  patterns  (which  were 
different  due  to  the  use  of  circular  polarization  which  required  ortho- 
gonally polarized  transmitter  and  receiver). 

The  radar  system  used  a Klystron  transmitter  and  low-noise  parametric 
receivers  which  could  be  cooled.  A coherent  linear  receiver  system  was 
used  with  matched  filters  on  both  the  sine  and  cosine  output  channels. 

The  output  from  each  filter  was  sampled  and  stored  in  real  time  in  a 
SDS  9300  computer.  The  data  were  square  law  detected  and  incoherently 
integrated  in  the  computer.  An  estimate  of  the  receiver  noise  power  was 
subtracted  from  the  detected  signal  plus  noise  to  obtain  a signal  power 
estimate.  Fifty  independent  samples  were  processed  resulting  in  a sta- 
tistical uncertainty  of  0.7  dB  rms  in  a back  s'-atter  cross  section  per 
unit  volume  estimate.  To  obtain  independent  pulses,  a sampling  rate  of 
20  jnilses  per  second  was  used  taking  2.5  seconds  to  obtain  an  observation. 

Radar  system  performance  was  continuous!}  monitored.  The  receiver 
system  was  continuously  monitored  using  a pulsed  noise  tube  of  known 
temperature.  The  average  transmitter  power  was  continuously  recorded. 
During  19'0  an  additional  calibration  pulse  was  used  to  provide  a second 
receiver  calibration.  The  two  independent  receiver  calibrations  differed 
by  less  than  -*0.2  dB. 

The  absolute  radar  system  calibration  was  made  by  combining  the 
values  for  each  factor  in  the  radar  equation.  The  maximum  error  in  the 
calibration  constant  was  estimated  at  ±1.4  dB  (at  3 standard  deviations). 
An  overall  calibration  check  was  made  by  tracking  orbiting  spheres  of 
known  cross  section.  The  measured  cross  sections  were  within  0.5  dB  of 
the  known  (true)  cross  sections.  Recent  observations  made  by  the  author 
in  1974  using  an  improved  computer  aided  satellite  tracking  system  and 
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TABLE  II 


Increased  to  80  dB  by  a manually  adjusted  variable  attenuator  in  the  transmission  line  between  the  antenna  and  receiver. 


the  Millstone  radar  produced  cross  sections  within  0.25  dB  of  the  true 
value.  Intercomparisons  between  Millstone  Hill  L-band  radar  measurements 
of  rain  and  simultaneous  measurements  with  the  independently  calibrated 
Avon-to-Westford  Bistatic  Radar  System  of  the  same  rain  volume  produced 
differences  of  less  than  1 dB  when  corrected  for  known  rain  drop  shape 
effects.  These  results  show  that  the  most  likely  absolute  calibration 
error  for  the  radar  system  is  in  the  ±0.4  - ±0.5  dB  range. 

4.2  Data  Processing 

4.2.1  Radar  Data 

The  radar  observations  were  recorded  on  magnetic  tape  for  post  test 

analysis.  The  data  were  processed  to  obtain  vertical  surface  distance, 

2 

height  (RHI)  maps  of  estimates.  These  values  were  further  integra- 
ted over  a 22.5  km  surface  distance  and  1 km  height  interval  at  a con- 

2 

stant  height  to  further  improve  the  signal  to  noise  ratio  of  a C pro- 

2 ” 

file  estimate.  The  calculated  minimum  single  pulse  value  detectable 
for  a non-fluctuation  target  at  a range  of  100  km  with  a unity  signal-to- 
noise  ratio  is  1 x 10  m as  indicated  in  Table  II.  Using  both 
noise  power  subtraction  in  the  process  of  signal  detection  and  the  addi- 
tional 15  range  cell  integration  implied  in  the  processing  to  obtain 
2 

a C profile  improved  the  signal-to-noise  ratio  by  14  dB  yielding  a 

^ 2 -18-2/3 

minimum  detectable  C value  of  4 x 10  m"  ' . In  practice,  the  mini- 

^ - 1 8 - 1 7 

mum  detectable  value  ranged  between  8 x 10  and  2 x 10  at  heights 

below  20  'm  depending  upon  the  system  configuration  and  the  extent  to 
which  the  system  was  bothered  by  man-made  interference. 

The  radar  system  was  bothered  by  sidelobe  detection  of  ground  clut- 
ter, aircraft,  and  heavy  rain.  To  circumvent  the  effect  of  ground  clut- 
ter, data  were  processed  only  for  ranges  in  excess  of  60  km.  CAT  detec- 
tion was  not  attempted  during  periods  of  heavy  rain;  radar  operations, 
however,  were  conducted  during  such  periods  because  the  investigation  of 
rain  scatter  was  the  principal  motivation  for  the  measurement  program 
(Crane,  1974).  The  effects  of  aircraft  and  ground  clutter  at  long  ranges 
during  periods  with  ducting  were  more  difficult  to  combat. 

Figures  1 and  2 depict  the  raw  radar  measurements  and  a computer 
generated  surface  distance,  height  map  (RHI).  In  figure  2,  the  intensity 
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Figure  2 Computer  generated  RHI  (note  the  correct  date  is  May  28,  19681 
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changes  represent  5 dB  steps  in  The  raw  radar  data  show  strong  echoes 

from  turbulence  in  the  planetary  boundary  layer  (height  less  than  3 km), 

2 -17 

an  intermediate  layer  of  weak  turbulence  'c-  10  , see  figure  4), 

solid  cirrus  overcast  in  the  7 to  9 km  height  region  and  4 turbulent 
layers  above.  The  layer  at  10  km  height  reveals  a billow-like  structure 
(at  the  resolution  of  this  radar)  with  strong  intensity  variations  at  an 
8 km  spacing.  Also  evident  in  the  data  are  enhanced  ground  clutter.  The 
strong,  constant  range  return  at  125  km  corresponds  to  the  location  of 
Cape  Code  as  viewed  from  the  Millstone  radar.  The  intermittent  arc  start- 
ing at  a surface  distance  of  65  km  is  typical  of  the  signature  of  an 
aircraft  in  a radar  sidelobe  when  tlie  range  to  the  aircraft  is  slowly 
changing  over  the  5 to  6 minute  time  period  reijuired  to  make  an  elevation 
scan . 

2 

I’rofiles  of  C^  were  constructed  for  3 separate  22.5  km  wide  surface 
distance  intervals,  77  to  99  km,  100  to  123  km,  and  124  to  146  km  respec- 
tively. The  profiles  obtained  for  the  radar  data  depicted  in  figures 
1 and  2 are  presented  in  figure  3.  The  three  profiles  are  combined  to 
make  one  best  estimate  profile  depicted  by  the  curve  labeled  1439  UT  on 
figure  4.  The  three  profiles  were  intercompared  to  remove  the  spurious 
effects  o1  aircraft  and  isolated  ground  clutter.  In  the  height  interval 

between  2(  and  30  km,  the  profile  for  123.6  to  146,1  km  shows  large  vari- 
2 

ations  in  C^  due  to  the  strong  ground  clutter.  Similar  variation,  but 
of  smaller  magnitude  is  seen  in  the  curve  for  100.2  to  122.7  km  and  no 
variation  is  seen  in  the  third  profile.  Spurious  effects  were  removed 
by  picking  the  minimum  value  from  the  three  curves  at  each  height. 

The  resultant  profile  is  displayed  in  figure  4.  Guidance  in  profile 
estimation  is  also  taken  from  the  RHI  display.  The  76.8  to  99.3  km  pro- 
file was  used  for  the  11  to  16  km  height  region  because  it  showed  equi- 
valent apparent  conlamination  to  the  minimum  value  profile  and  better 
reproduced  the  structure  evident  in  the  65  to  80  km  surface  distance 
interval  at  those  heights.  As  seen  from  the  profile  data  in  figure  3, 
the  differences  between  the  two  possible  profiles  in  the  11  to  16  km 
height  region  are  less  than  1 dB. 

The  apparent  spread  in  the  width  of  the  cirrus  cloud  layer  in  the 
computer  generated  RHI  is  due  in  part  to  the  method  used  to  construct 
the  RHI.  Radar  data  from  each  range,  elevation  interval  was  assumi;d  to 
uniformly  fill  the  radar  resolution  volume.  The  radar  beamwidth  is  0.7° 
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spanning  a 1.2  km  height  interval  at  100  km  range.  When  the  layer  is 
thinner  than  the  beamwidth,  the  amplitude  of  the  layer  changes  with  dif- 
ference to  the  profile.  When  the  layer  is  thicker  than  a beamwidth,  as 
in  this  case,  the  magnitude  is  constant  but  the  width  appears  to  increase 
with  distance  to  the  layer.  When  the  profiles  are  nearly  identical  as 
in  the  case  considered  above  for  the  11  to  16  km  height  region,  the  pro- 
file for  the  smallest  surface  distance  is  used  because  of  its  better 
height  resolution. 

2 

The  computer  processing  used  to  produce  the  profiles  presented 

in  figure  3 was  accomplished  at  Lincoln  Laboratory  between  1970  and  1972. 

The  data  tapes  used  to  generate  the  plots  were  released  during  the  1974- 

1975  time  period  when  support  for  further  processing  was  not  available. 

For  this  contract  further  processing  was  accomplished  manually  using  the 

2 

comp  Iter  generated  profiles.  A computer  profile  was  generated  for 

each  RHI  scan.  The  median  profile  for  5 successive  RHI  scans  was  obtained 

2 

as  a best  estimate  of  the  mean  profile  for  a 20  to  30  minute  period. 

The  median  of  these  profiles  was  then  constructed  as  a best  estimate  for 

the  mean  of  an  observation  set.  The  median  of  the  medians  is  used  to 
2 

characterize  the  C profile  for  one  set  of  observations.  It  is  used  to 
n 

provide  the  representation  of  the  time  average  of  the  profile  required 
for  use  in  estimating  the  eddy  parameters  as  described  in  section  2.  The 
median  and  two  sample  profiles  taken  6 minutes  apart  are  shown  in  figure 
4.  Although  the  data  for  1439  and  1445  UT  are  obtained  from  spatial 
locations  nearly  140  km  apart,  the  profiles  in  the  region  of  turbulence 
above  a height  of  10  km  are  nearly  identical.  The  height  of  the  cirrus 
cloud  layer  does  change  over  that  spatial  distance. 

The  median  for  the  observation  set  is  roughly  a factor  of  4 lower 
at  12-14  km  height.  The  receiver  noise  line  displayed  on  the  figure 
represents  the  expected  median  profile  if  receiver  noise  alone  were  present. 
The  variation  of  the  median  curve  about  this  line  at  heights  above  about 
17  km  may  be  used  as  an  indication  of  the  statistical  uncertainty  in  the 
median  estimate.  Where  possible,  median  estimates  were  obtained  only 
for  observations  in  a narrow  azimuth  sector  of  less  than  40“  to  remove 
any  possible  large  scale  spatial  changes.  For  the  data  on  figure  4, 
neither  single  observation  profile  was  used  in  the  construction  of  the 
median.  This  selection  process  did  not  appear  to  change  the  resultant 
median  suggesting  that  the  steady  state,  horizontally  homogeneous  assump- 
tions are  valiil  on  the  time  scale  used  for  analysis. 
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The  median  profiles  for  each  observation  set  are  given  in  Appendix 


B. 


4.2.2  Radiosonde  Data 

Standard  radiosonde  observations  were  specially  made  by  the  Air 
Force  Cambridge  Research  Laboratory  (now  Air  Force  Geophsyics  Laboratory) 
in  support  of  the  radar  measurement  program.  These  data  were  computer 
processed  at  MIT  Lincoln  Laboratory  using  raw  data  samples  at  30  second 
intervals  in  addition  to  the  significant  level  values.  On  several  ocas- 
sions,  Loran  C wind  sensors  were  flown  on  the  balloons  to  provide  im- 
proved wind  velocity  profiles.  The  data  were  processed  to  provide  poten- 
tial temperature,  refractive  index,  and  windspeed  and  direction  profiles. 
Sample  profiles  for  the  data  presented  in  figures  1 to  4 are  presented 
in  figures  5 to  9.  Two  sets  of  soundings  were  available  on  many  of  the 
days.  Although  the  water  vapor  contribution  to  the  refractive  index  was 
not  explicitly  treated  in  Section  2,  the  data  were  available  for  analysis 
as  shown  in  figure  8.  Figure  9 presents  the  radio  reflectivity  profile 
where  refractivity  is  (n-1)  x 10^  (see  equation  (7)). 

The  radiosonde  data  were  processed  in  1 km  height  intervals.  The 
layer  thickness  h was  taken  to  be  1 km  to  approximate  the  vertical  resolu- 
tion of  the  radar  system  at  the  70  to  140  km  ranges  and  the  low  elevation 
angles  used  for  the  measurements.  The  data  required  to  calculate  c^, 
k'lj,  and  d were  selected  from  the  radiosonde  data  computer  output  (print 

out)  and  used  to  manually  calculate  the  parameters  (using  a desk  calcu- 

_ 2 

lator).  Intermediate  profiles  of  were  generated  and  are  separately 

displayed  as  shown  in  figure  10.  Profiles  of  and  e were  then  cal- 

n in  2 

culated  as  displayed  in  figure  11.  A complete  set  of  C.j,  , d,  and 
profiles  are  given  in  Appendix  A. 
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Figure  5 Potential  temperature  profile,  May  28,  1968 
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5.  APPLICATION  OF  THE  TURBULENCE  MODEL 


5.1  Comparison  with  U-2  Pilot  Reports 

The  HICAT  U-2  aircraft  was  stationed  at  Hanscom  AFB  during  May 
through  July  1968,  for  equipment  installation.  During  this  interval, 
arrangements  were  made  for  7 U-2  missions  on  days  the  radar  was  being 
used  for  CAT  measurements.  The  HICAT  instrumentation  was  not  available 
for  these  flights  but  the  pilots  had  been  associated  with  the  HICAT  pro- 
gram and  had  some  experience  in  "seat  of  the  pants"  detection  of  CAT. 

The  aircraft  was  flown  on  radial  paths  relative  to  the  radar  and 
used  to  sound  the  atmosphere  between  10  and  20  km  height.  The  pilot 
subjectively  recorded  the  intensity  of  turbulence  and  the  height  inter- 
val over  which  it  was  observed.  Pilot  reports  ranged  from  extremely 
light  (EL)  to  moderate  (M) . The  radar  derived  estimates  for  and  the 
pilot  reports  are  depicted  in  figures  12-16.  The  radar  derived  eddy 
dissipation  rate  estimates  include  both  the  mean  estimated  layer  values 
and  active  layer  values.  In  addition,  the  highest  observed  single  scan 
values  for  the  azimuth  and  time  interval  of  the  aircraft  measurements 
are  displayed. 

A calibration  of  the  pilot  is  required  to  make  a comparison  with 
the  radai  data.  The  pilots  employed  the  same  intensity  classification 
procedure  they  used  during  the  HICAT  program.  During  HICAT,  instrumenta- 
tion was  used  only  for  turbulence  catalogued  as  light  or  more  intense. 

The  pilot  reports  were  recorded  for  each  of  the  data  samples  obtained 
during  the  HICAT  program.  For  example,  spectra  obtained  from  two  passes 
through  a region  of  turbulence  are  displayed  in  figure  17  (obtained  from 
Crooks,  it  al,  1967,  Vol.III).  These  spectra  are  for  turbulence  reported  as 
light  by  the  pilot.  The  reported  wavenumbers  are  for  large  scale  aniso- 
tropic turbulence  and  do  not  reveal  the  existence  of  an  isotropic  iner- 
tial subrange.  In  the  isotropic  inertial  subrange  all  the  spectra  as 
scaled  for  display  on  this  figure  should  lie  one  on  top  of  the  other. 

The  equii'alencc  of  the  spectra  for  the  two  passes  again  indicates  the 
utility  I f the  steady  state,  horizontal  homogeneity  assumptions  for 
stratospheric  turbulence.  The  short  straight  line  superimposed  on  the 
figure  is  an  attempt  to  estimate  the  behavior  of  the  spectra  in  the  iso- 
tropic inertial  subrange.  The  eddy  dissipation  rate  value  for  these 
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Figure  15  Hadar  and  Aircraft  comparison;  e profiles,  July  30,  1968 
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Figure  16  Radar  and  Aircraft  comparison;  e profiles,  July  31,  1968 
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Figure  17  One-dimensional  velocity  spectra  obtained  from  Hli.AT  (Crook 
et  al,  1967) 
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spectra  was  calculated  using  the  location  of  this  approximate  line  on 

the  plot  (using  a power  spectral  density  value  at  a given  wavenumber) . 
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A value  of  2 x 10  m /s  was  obtained  as  reported  on  the  figure.  This 
process  was  repeated  for  a number  of  spectra  catalogued  as  corresponding 
to  light  or  moderate  turbulence.  The  resultant  e values  were  used  to 
position  the  L (light)  and  M (moderate)  values  on  figures  12-16.  The 
calibration  of  EL  (extremely  light)  through  VL  (very  light)  could  not 
be  made  using  reported  spectra  since  they  correspond  to  turbulence  of 
insufficient  intensity  to  produce  reliable  estimates  of  the  spectra. 

These  levels  were  approximated  using  information  reported  by  Crooks  et 
al  (1967)  and  supplied  by  the  pilots. 

The  U-2  reports  are  for  single  passes  through  active  regions  of 
turbulence  exceeding  the  intensity  - layer  thickness  threshold  needed  to  ! 

trigger  a pilot  response.  Since  the  maximum  to  median  ratios  for 

typically  varied  by  less  than  a factor  of  4,  it  is  expected  that  the  | 

maximum  to  median  e values  should  vary  by  less  than  a factor  of  8.  i 

Agreement  between  the  pilot  reports  (uncertain  by  a factor  of  4 perhaps)  j 

and  the  radar  derived  estimates  is  obtained  when  the  pilot  reports  are  i 

within  a factor  of  10  of  c and  a factor  of  4 of  e for  the  single  scan.  j 

The  U-2  reports  should  indicate  more  intense  turbulence  than  the  radar. 

This  behavior  is  apparent  for  all  the  U-2/radar  comparisons  with  the 
exception  of  obseirvations  made  below  11  km  on  29  May  and  5 June  (figs. 

12  and  13).  The  radar  data  for  29  May  are  probably  contaminated  by  cir- 
rus clouds.  This  possibility  also  exists  for  5 June  although  no  visual  J 

reports  of  cirrus  were  received.  The  results  of  this  comparison  indi-  i 

cate  that  within  the  uncertainties  associated  with  the  use  of  pilot 
reports,  excellent  agreement  was  obtained. 

The  U-2/radar  comparisons  are  the  most  extensive  collection  of  data 
available  for  the  validation  of  the  model,  equations  (26)  and  (27).  The 
only  other  available  data  useful  for  model  verification  were  reported  by 
Metcalf  and  Atlas  (1973).  They  reported  maximum  radar  cross  section 
values  obtained  using  a radar  resolution  volume  smaller  than  the  region 
of  active  turlmlence.  They  also  reported  aircraft  observations  made 
through  the  turbulent  region.  Their  turbulent  energy  balance  however, 
was  calculated  based  upon  data  from  a region  thicker  than  the  peak  of 
the  turbulent  layer.  Based  upon  their  data  and  equation  (27),  the  eddy 
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dissipation  rate  at  the  peak  of  their  layer  should  be  in  the  10-40  cm  /s 

2 3 

range.  They  report  estimated  eddy  dissipation  rates  in  the  4-26  cm  /s 
range,  the  value  depending  upon  the  procedure  they  used  for  estimation. 

Considering  the  uncertainties  in  their  data  particularly  with  respect  to 

the  relative  effects  of  the  potential  temperature  and  absolute  humidity 

gradients  and  the  discrepancy  associated  with  using  the  maximum  radar 

return  with  a 1 meter  (height)  resolution  volume  to  compare  against  air-  ■ 

craft  data  with  a 20  m (height)  resolution  volume,  excellent  agreement  is 

apparent. 

5.2  Eddy  Dissipation  Rate  Estimates 

Median  profiles  of  eddy  dissipation  rate  were  obtained  for  each 
observation  set.  The  individual  profiles  are  presented  in  Appendix  A. 

Composite  profiles  organized  by  observing  period  are  presented  in  fig- 
ures 18  through  24.  These  figures  show  rather  large  variations  from  one 
observation  set  to  another  at  heights  below  about  12  km  and  little  vari- 
ation above.  In  the  region  between  12  and  16  km  most  of  the  observed 
median  values  are  within  an  order  of  magnitude  range.  These  results 
indicate  that  although  a rather  large  variation  in  turbulent  energy  dissi- 
pation may  occur  in  the  lower  troposphere,  the  average  dissipation  rate 
within  the  barrier  between  the  troposphere  and  stratosphere  is  remarkably 
consistant  from  day-to-day. 

Data  that  were  obviously  contaminated  by  scattering  from  cirrus 
clouds  were  not  included  in  these  figures.  Some  contamination  still  may 
be  present  at  heights  below  12  km.  The  minimum  on  average  at  the  barrier  ; 

is  real.  The  median  of  ail  the  medians  was  used  to  obtain  the  values  j 

given  in  Table  I.  To  reduce  possible  biasing  affects  caused  by  cirrus  j 

clouds  at  heights  below  12  km  and  by  radar  sensitivity  at  heights  above  I 

16  km,  values  at  these  heights  were  not  entered  in  Table  I.  | 

5.3  Eddy  Diffusion  Coefficient  Estimates 

The  median  eddy  diffusion  parameters  are  presented  in  figures  25 
through  31 . The  general  shapes  of  the  curves  are  identical  to  the  eddy 
dissipation  rate  curves.  These  data  demonstrate  the  existence  of  a 
barrier  in  the  lower  stratosphere  to  vertical  transport  by  small  scale 
turbulent  processes.  The  increase  of  above  the  barrier  is  in  agreement 
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Figure  18  Eddy  dissipation  rate  profiles,  winter  1968 
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Figure  19 


Eddy  dissipation  rate  profiles,  spring  1968 
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Figure  20  Eddy  dissipation  rate  profiles,  summer  1968 
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Figure  21  Eddy  dissipation  rate  profiles,  summer  1970 
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Figure  22  Eddy  dissipation  rate  profiles,  summer  1970 
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Figure  24  Eddy  dissipation  rate  profiles,  fall  1970 
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with  the  shape  of  the  eddy  diffusion  curves  used  in  large  scale  trans- 
port models  (see  Reiter,  1975),  Radar  observations  of  vertical  transport 
coefficients  in  the  mesosphere  indicate  that  the  values  continue  to 
increase  in  height  eventually  merging  with  the  values  assumed  by  larger 
scale  modelers  (Cunnold,  1975)  indicating  that  the  small  scale  turbulent 
processes  are  the  most  important  processes  at  the  higher  altitudes  well 
above  the  lower  stratosphere  barrier. 

5.4  Active  Layer  Thickness  Estimates 

The  model  developed  for  the  estimation  of  and  using  the  radar 

data  required  an  estimate  of  d since  the  dominant  region  of  active  tur- 
bulence is  thinner  than  the  resolution  volume  of  the  radar.  If  this 
analysis  were  not  performed,  the  and  values  would  be  underestimated 
since  the  turbulence  did  not  "fill  the  beam".  The  active  layer  thick- 
ness estimates  were  obtained  directly  from  radiosonde  data  assuming  only 
that  the  atmosphere  was  everywhere  turbulent  but  with  no  knowledge  of 
the  intensity  of  the  turbulence.  Only  the  aircraft  or  radar  can  provide 
estimates  of  the  intensity  of  turbulence. 

The  median  layer  thickness  for  all  heights  in  the  12  through  16  km 

region  was  estimated  using  the  radiosonde  data  (all  e values),  all  the 

^ 2 

available  Lf-2  observations  above  2 x 10  m /s  or  all  pilot  reports), 

and  for  all  U-2  observations  that  would  have  been  processed  using  the 

-4  2 3 

Project  HICAT  criterion  (e  above  2 x 10  m /sec  or  all  light  and  mod- 

mm  ° 

erate  turbulence).  The  results  are  shown  in  Table  III  together  with  the 
number  of  samples  used.  These  data  show  that  the  median  turbulent  layer 
size  is  thin  and  the  stronger  layers  are  very  thin.  It  is  interesting 
to  note  that  these  observations  are  at  variance  with  the  results  reported 
by  Crooks  et  al  (1967)  which  showed  a median  layer  thickness  of  750  m 
for  13  samples,  a factor  of  15  times  thicker  than  our  results.  Our  data 
were  obtained  over  flat  terrain  but  no  indication  was  given  about  the 
conditions  under  which  the  HICAT  data  were  obtained. 

5.5  Temporal  Variation  of  Turbulent  Layer  Intensity 

The  raw  elevation  scan  data  - the  RHI's  such  as  Figure  1 - always 
depict  CAT  as  occurring  in  thin  layers.  The  weaker  layers  are  sometimes 
evident  only  in  a single  RHI  observation;  the  stronger  layers  may  persist 
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TABLE  III 

ACTIVE  LAYER  THICKNESS  ESTIMATES 


d{m) 

# of  samples 

All  !-ayers  from  Radiosonde  Data 

200 

236 

All  11-2  Observations 

100 

49 

All  U-2  Data  in  the  Light  to 
Moderate  Category 

50 

12 

‘ for  one-to-two  hours.  The  layer  at  the  tropopause  may  persist  through- 

out an  entire  measurement  period.  The  temporal  variation  of  the  layer 

f structure  for  28  May  1968  is  revealed  in  Figure  4 and  Figure  32.  The 

2 2 
median  profile  in  Figure  4 shows  a general  decrease  in  with  height  | 

! above  the  tropopause  (''-12  km).  The  two  "instantaneous"  profiles  - RHI's  ; 

[ taken  6 minutes  apart  - show  a definite  layer  structure.  The  time  history 

f of  the  layer  structure  is  revealed  in  Figure  32.  The  data  used  to  gener-  j 

1 ate  this  figure  were  obtained  from  the  raw  RHI  photographs  and  the  layer 

i intensities  coded  on  the  figure  are  only  approximate. 

The  time  history  reveals  a persistant  layer  at  the  tropopause.  The  J 

j height  of  this  layer  varies  over  the  observing  period  producing  a smear- 
ing of  the  structure  as  revealed  by  the  median  profile.  The  layers  above  |‘ 

the  tropopause  persist  for  one-to-two  hours.  The  layers  in  this  region  !: 

2 ■ ' 

seem  to  occur  at  different  heights.  The  general  decrease  in  with 

height  depicted  by  the  median  profile  reflects  the  variation  in  height  ; 

of  the  persistant  layer  at  the  tropopause  and  the  random  occurrence  of 

the  layers  at  higher  altitudes. 

2 

The  median  profiles  depicted  in  the  figures  in  Appendix  B reveal 
the  same  behavior.  Median  profiles  for  one-to-two  hour  time  durations 
sometimes  show  a layer  structure  as  evident  at  a height  of  19  km  on  the 
13  February  '68  profile  but  at  other  times  show  little  structure  indica- 
tive of  weak,  randomly  occuring  layers  (see  the  19  February  '68  profile) . 
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6.  CONCLUSIONS 

A model  was  developed  for  the  combination  of  radar  and  radiosonde 
data.  The  basic  premise  used  for  the  development  of  the  model  was  that 
the  atmosphere  and  all  its  sublayers  is  always  turbulent.  The  intensity 
of  the  turbulence  may  vary  by  orders  of  magnitude  from  layer  to  layer. 

The  mean  profiles  provided  by  radiosonde  measurements  cannot  be  used  to 
estimate  the  intensity  of  turbulence  - the  eddy  dissipation  rate  - but 
can  only  be  used  to  characterize  the  stability  of  the  layers.  Additional 
data  must  be  provided  to  estimate  the  intensity  of  turbulence.  These 
data  may  be  obtained  from  well  calibrated  radars  or  from  instrumented 
aircraft.  In  the  analysis  performed  by  VanZandt  et  al  (1977)  the  re- 
quired additional  data  were  assumed.  In  the  analysis  performed  by 
Rosenberg  and  Dewan  (1974),  the  additional  data  were  also  assumed,  this 
time  by  picking  a value  for  the  time  step. 

The  model  developed  for  the  analysis  of  radar  and  radiosonde  data 
was  shown  to  be  consistant  with  all  the  available  measurements.  The 
results  match  closely  results  obtained  by  Lilly  et  al  (1974)  for  similar 
conditions  of  underlying  terrain.  The  data  show  the  existence  of  a 
barrier  to  vertical  transport  by  small  scale  turbulence  between  the 
troposphere  and  stratosphere.  The  small  scale  turbulent  transport  pro- 
cess over  flat  land  at  a midlatitude  location  is  an  order  of  magnitude 
smaller  than  required  to  explain  the  observed  distributions  of  trace 
constituents  using  one-dimensional  models.  This  discrepancy  may  be 
caused  either  by  the  neglect  of  larger  scale  transport  processes  or  by 
strong  geographical  variations  in  the  intensity  of  turbulence  within  the 
barrier.  The  results  presented  above  show  that  seasonal  variations  can- 
not explain  the  discrepancies. 
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Figure  10  C , C and  active  layer  thickness  profiles,  May  28,  1968 
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